Top-Quark Mass Measurement in the Dilepton Channel Using in situ Jet Energy Scale 

Calibration 
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We employ a top-quark mass measurement technique in the dilepton channel with in situ jet 
energy scale calibration. Three variables having different jet energy scale dependences are used 
simultaneously to extract not only the top-quark mass but also the energy scale of the jet from 
a single likelihood fit. Monte Carlo studies with events corresponding to an integrated luminosity 
of 5 fb _1 proton-proton collisions at the Large Hadron Collider y/s — 7 TeV are performed. Our 
analysis suggests that the overall jet energy scale uncertainty can be significantly reduced and the 
top-quark mass can be determined with a precision of less than 1 GeV/c 2 , including jet energy scale 
uncertainty, at the Large Hadron Collider. 
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In the standard model, the top quark (i) is the heaviest 
known elementary particle [lj. The heavy mass signifi- 
cantly affects the electroweak radiative correction that 
relates the top-quark mass and the W boson mass to 
the Higgs boson mass |2j. In addition, the heavy mass 
may have implications for new physics theories, includ- 
ing the minimal supersymmetric standard model and 
technicolor-like models. Because of its importance, top- 
quark mass (Mt p) measurements have been performed 
using various methods in different decay channels. The 
precision of the M top measurement already surpasses 
0.5% at the Tevatron 0. 

It is important to measure Mt op using different tech- 
niques and independent data samples in different decay 
channels. Significant differences in the measurements of 
Mtop using different decay channels may indicate contri- 
butions from new physics beyond the standard model [1] . 
The dilepton decay channel is particularly interesting be- 
cause the signature of this channel can be mimicked with 
supersymmetric partner stop pairs Q as well as charged 
Higgs boson signals @. However, the M top precisions of 
the dilepton channel measurements have been lim- 
ited because the branching fraction is much smaller and 
systematic uncertainty from the jet energy scale (JES) is 
much larger than those of the lepton+jets channel mea- 
surements [Io|-[l3|]. At the Tevatron, a single measure- 
ment already achieved 1.2 GeV/c 2 precision in the lep- 
ton+jets channel [l(| ■ However, the total uncertainty in 
the dilepton channel in a standalone measurement was 
3.0 GeV/c 2 8], with the uncertainty from JES being 
dominant. In the Large Hadron Collider (LHC), the 
cross section of ti pair production is approximately 20 
times larger than that in the Tevatron. Therefore, the 
small branching fraction corresponding to data obtained 
in LHC experiments up to 2011 (more than 5 fb _1 inte- 



tainty, (i^'fj GeV/c 2 ), which is dominated by the overall 
JES uncertainty (2.0 GeV/c 2 ). The JES systematic un- 
certainty cannot generally be reduced by using a larger 
data sample. 

A similar issue was encountered in the lepton+jets 
channel. However, in situ JES calibration using hadronic 
decaying W bosons [HI resolved this issue. In this 
method, the overall JES uncertainty was absorbed in the 
statistical uncertainty; therefore, a larger data sample 
could reduce not only the pure statistical uncertainty but 
also the overall JES systematic uncertainty. However, in 
the dilepton channel, no single variable can be used to 
calibrate the JES uncertainty in situ for the M top mea- 
surement. A recent measurement performed by the DO 
Collaboration employed a new technique for JES cali- 
bration in this channel [HI]. They used the result of the 
JES measurement in the lepton+jets channel. In this 
way, they significantly reduced overall JES systematic 
uncertainty in the dilepton channel. However, it is not a 
standalone measurement and the use of the lepton+jets 
channel result introduced an additional systematic un- 
certainty. 

In this Letter, we propose a novel technique for the 
Mtop measurement in the dilepton channel, in which we 
perform the in situ JES calibration using three variables. 
The selected variables, which have already been used for 
the M top measurements, have different dependences on 
M top and JES, so we can perform simultaneous measure- 
ment of both M t0 p and JES. This is the first time the 
M t0 p measurement has been performed using in situ JES 
calibration in the dilepton channel even though it is a 
Monte Carlo (MC) simulated experiment. For a realistic 
estimation of the precision, we consider the environment 
of the LHC experiment with 5 fb _1 pp collisions. 

There arc different definitions of quark masses in the 



grated luminosity) is not an issue. As an example, the theoretical framework [16j. M top of the MS renormal- 



recent M top measurement in the dilepton channel per- 
formed by the CMS Collaboration had 1.2 GeV/c 2 sta- 
tistical uncertainty using 2.2 fb _1 data @. However, 
this measurement returned a larger systematic uncer- 



ization scheme, M^s, differs from the pole mass, M p i oc , 
by about 10 GeV/c [la, The direct M top measure- 
ments are all calibrated with MC simulations; the mea- 
sured quantities therefore correspond to the M t0 p scheme 
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used in the MC simulations, Mmc- Even though it is usu- 
ally assumed that M M c is the same as M pole Q, there 
are a number of theoretical questions in relating M po i e 
to Mmc [3 13 • An accurate relation is still under the- 
oretical investigation [l9j]. As in all other direct M top 
measurements, we measure Mmc in this Letter. 

We generate simulated ti samples using the leading 
order (LP) MC generator MADGRAPH/madevent pack- 
age [20]. In the madgraph/madevent generation, we 
vary the parameter M top between 160 and 190 GeV/c 2 in 
steps of 2 GeV/c 2 . A total of 500 000 ti pair events in the 
dilepton final state are produced for each sample. Show- 
ering and hadronization are performed by pythia plj ]. 
To take into account the detector effect and perform 
event reconstruction under realistic conditions, we use 
the fast detector simulation package delphes 22]. The 
conditions used for the simulated detector are the ones 
generally associated with the LHC detector. We assume 
the coverage of the tracker to be within \rf\ — 2.5 with 
100% efficiency and that of the calorimeter to be within 
\r)\ = 3 with tower segment A77 ~ 0.1 and A(f> ~ 0.1. The 
resolution of the electromagnetic calorimeter (EM) and 
hadronic calorimeter (HA) are parameterized by 
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All physics objects such as leptons, jets, and missing 
transverse energy are reconstructed in the fast simula- 
tion. Jets originating from b quarks are tagged using a 
reliable 6-tagging algorithm with an efficiency of approx- 
imately 40%. 

In the dilepton decay channel, the production of a 
ti pair is followed by the decay of each top quark to 
a W boson and a b quark, where both W bosons de- 
cay to charged leptons (electron or muon) and neutri- 
nos (ti — > ll'vv'bb). Events in this channel thus con- 
tain two leptons, two b quark jets, and two undetected 
neutrinos. To select the candidate events of ti dilepton 
topology, we require two oppositely charged lepton can- 
didates with pt > 20 GeV/c. We also require a miss- 
ing transverse energy exceeding 25 GeV and at least two 
tagged jets with Et > 30 GeV. The expected signal and 
background events are 5961 ± 545 and 320 ± 97, respec- 
tively, taken from Ref. @ with scaling according to the 
respective integrated luminosity for the 5 fb _1 data. The 
expected background contribution is approximately 5%. 
Because of its small contribution, we do not include back- 
ground events in our further analysis, for simplicity. 

We use three reconstructed variables as observables of 
M top as well as the nuisance parameter Ajes- Ajes is 
the relative energy scale of the jet with respect to the 
nominal JES as used in the lepton+jets channel mea- 
surements [Icil - [l3j . All three variables have already been 
used in the M top measurements. 



using the neutrino weighting algorithm [23j: 
is still widely used in the dilepton channel 0-0] as 
a template method. A study by the CDF Collab- 
oration indicates that m^ w has the most-precise 
statistical uncertainty among the variables used in 
its study [11]. 

mT2 - The collider variable related to the missing 
transverse mass of the system of two missing par- 
ticles: tot2 was initially developed for new physics 
particles as a variable sensitive to the mass of new 
particles in pair production (25[. Because of the 
two missing particles (neutrinos) in the dilepton 
channel, this variable was suggested for the M top 
measurement [2q . The CDF Collaboration per- 
formed the M top measurement using 777x2 and im- 
proved the precision with the simultaneous usage of 
mf w and 777x2 [2J]- They also extensively studied 



the statistical as well as systematical uncertainties 
in the different variables. The expected statisti- 
cal uncertainty obtained using m^ w had the small- 
est value; however, the systematic uncertainty was 
larger than that in the measurement performed us- 
ing 771T2 ■ This difference was caused by different 
JES systematics due to the different dependence of 
each variable on M top and JES. 

o lepton rTHi r j. 1 j. lepton 

3. jj t - Ihe average px of two leptons: p T 
has been suggested as a good variable for the M top 
measurement with large statistics in the LHC 27]. 
jjy Pton , in general, does not depend on JES, and 
it negligibly contributed to the systematic uncer- 
tainty from JES. The CDF Collaboration has per- 
formed measurements using this variable in the lep- 
ton+jets [HI and the dilepton [29J channels, ob- 
taining negligible JES systematic uncertainty but 
relatively large statistical uncertainty. This is basi- 
cally caused by the low sensitivity of M top and the 
insensitivity of JES to the lepton px variable. 

The three variables discussed above have very different 
A/top and JES dependences. Therefore, it is possible to 
extract both Aft op and JES information together if we use 
the three variables simultaneously. To simulate and ex- 
tract JES information from the MC pseudoexperiments, 
we vary the scale of the jet energy relative to nominal 
JES (Aj ES ) in the MC simulations from -15% to 15% 
in 1.5% steps. Even though the overall JES uncertainties 
depend on the jet 77 and px, we use the overall variation 
for simplicity in this study. The overall jet energy uncer- 
tainty in the CMS experiment [3(| is expected to be from 
2% to 3%. The plots of each variable are shown in Fig. [T] 
with different M top as well as Ajes- mf w and mx2 de- 
pend on both M top and Ajes- However, it is clear that 
^tepton ^ oes nQt ^ e p en( j on Ajes- Therefore, using three 
variables simultaneously in a single likelihood function, 
we can extract Ajes information in situ as a nuisance 
parameter of the M top measurement. 
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FIG. 1: Plots of the three variables for different M top (a, b, c) and Ajes (d, e, f). 



The technique of using three variables in a single like- 
lihood function was developed by the CDF Collabora- 
tion and has been used in lepton+jcts channel measure- 
ments. We follow the procedure discussed in Ref. Q. 
We estimate the probability density functions (PDFs) of 
signals, P sig (mf w ,m T2 ,^ pton ; M top , A, ms ), using ker- 
nel density estimation [31]. For the discrete values of 
M top and Ajes, we estimate the PDFs for the observ- 
ables. We smooth and interpolate the MC distributions 
to find PDFs for arbitrary values of M top and Ajes us- 
ing the local polynomial smoothing method [HJ . We then 
build the unbinned maximum likelihood [33[ for N events 
(and n s expected signal events) with the Poisson fluctu- 
ation: 

-n s N 

r c "« nJV n / NW lcpton 7i r a \ 

£ = — j-p — T, i=1 P aig (m t ,m T2 ,PT I M top , Ajes)- 

The quantity P slg (mf w , m T 2,PT Pt ° n ; M top , A JE s) de- 
notes the signal PDFs as determined by kernel density es- 
timation and local polynomial smoothing as a function of 
M top and Ajes- We minimize the negative logarithm of 
the likelihood using minuit [34j with respect to all three 
parameters (n s , M top , and Aj E s). The uncertainty on 
-Mop or Ajes is found by searching for the points where 
the negative logarithm of the likelihood minimized with 
respect to all other parameters deviates by 1/2 from the 
minimum. The uncertainty on the Af top measurement 
obtained in this way includes the statistical uncertainty 
as well as the systematic uncertainty of the overall JES 
owing to the allowed variation of Ajes ■ We also perform 
the likelihood fit without varying Aj E s (M top -only fit) 
and set it to zero. The M top -only fit measurement allows 
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FIG. 2: Negative log-likelihood contours of a single pseudo- 
experiment using a M top = 173 GeV/c 2 and Ajes = 0.0% 
sample. The minimum is indicated by the "x" and corre- 
sponds to the most probable top-quark mass and Ajes- The 
contours are drawn at 68% confidence level and 95% confi- 
dence level, corresponding to la and 2a uncertainty of the 
measurement, respectively. 



us to compare the performance of the new technique with 
the in situ JES calibration (2D fit). 

We test the likelihood procedure using MC pseudo- 
experiments. We construct pseudodata from a certain 
value of M t0 p and Ajes- We select the number of signal 
events from a Poisson distribution with a mean equal to 
the expected number of signal events, 5961 ± 545. We 
perform the maximum likelihood fit described in the pre- 
vious section. In principle, the likelihood fit, on average, 
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returns the value of the top quark mass used to gener- 
ate the pseudoexperiments. Figure [5] shows an exam- 
ple of a likelihood fit contour from a single pseudoex- 
periment of the 2D fit. We use a M top = 173 GeV/c 2 
and Ajes = 0% sample for this experiment and obtain 
Mtop = 172.89 ±0.57 GeV/c 2 and A JES = 0.22 ±0.64%. 
We verify the M top -only fit using the same pseudodata 
and obtain M top = 173.12 ± 0.37 GeV/c 2 . We perform 
this experiment 3 000 times for seven different M top val- 
ues ranging from 168 to 178 GeV/c 2 . Figures [3] (a) 
and (b) show the average residual (deviation from input 
Mtop) and the width of the pull (the ratio of the residual 
to the uncertainty reported by MINUIt), respectively, in 
the 2D fit using samples of Ajes = 0.0% without correc- 
tions. The small positive bias, 0.43 GeV/c 2 , is corrected 
and the uncertainty is correspondingly increased by 5% 
to correct for the width of the pull distribution. The 
residual and pull width for the Ajes parameter are also 
investigated using the same pseudoexperiments. We ob- 
tain a small negative bias of 0.33% and a pull width of 
1.08. We apply suitable corrections for the nuisance pa- 
rameter Ajes- These corrections enable us to test the 
M top measurement with different Ajes parameters. We 
vary the input Ajes from —3.0% to 3.0%, corresponding 
to approximately 1 to 1.5 times the JES uncertainty in 
the LHC experiment [3(| • Figure[3] (c) shows the residual 
distribution for various Ajes samples. There is no sig- 
nificant effect of Ajes on the mass residual. Therefore, 
we do not apply a correction of M top for the Ajes pa- 
rameter. The same procedures for the M top -only fit are 
performed, and we find no bias. The width of the pull 
distribution is also consistent with unity. 

With the corrections of the residual and the width of 
the pull, we can obtain the expected uncertainty. In 
Fig. 21 the expected statistical uncertainty for the M top - 
only fit (a) and the 2D fit (b) have been plotted. Because 
the 2D fit absorbs the overall uncertainty from JES, it has 
a larger statistical uncertainty than the M top -only fit. If 
we choose M t0 o = 173 GeV/c 2 , close to the world av- 
erage of M top 3, the 2D fit has an expected statistical 
uncertainty of 0.60 ± 0.03 GeV/c 2 , whereas that of the 
Mto P -only fit is 0.36 ±0.02 GeV/c 2 . However, the Mtop- 
only fit may have a much larger JES systematic uncer- 



tainty. To obtain the systematic uncertainty of JES, we 
check the mass residual as a function of Ajes, as shown 
in Fig. I (c) (Mto P -only fit) and (d) (2D fit). If we con- 
sider an optimistic 2% overall JES uncertainty [3(|, the 
Mtop-only fit gives an expected JES systematic uncer- 
tainty of 1.68 GeV/c 2 , whereas that in the case of the 
2D fit is 0.08 GeV/c 2 . If we consider a realistic JES 
systematic uncertainty, which will be 2%-3% depend- 
ing on jet px and n, the M top -only fit may have a JES 
systematic uncertainty slightly larger than 1.68 GeV/c 2 , 
which is consistent with the recent CMS measurement of 
2.0 GeV/c 2 without in situ JES calibration. However, 
the JES systematic uncertainty for the 2D fit will still 
be very small. We calculate the expected uncertainty 
from statistics and overall JES systematics together with 
a quadrature sum. The M top -only fit gives an expected 
uncertainty of 1.72 GeV/c 2 , whereas that in the case of 
the 2D fit is much smaller, 0.61 GeV/c 2 . The expected 
statistical uncertainty of Ajes m the 2D fit is 0.61±0.05% 
if we use M top = 173 GeV/c 2 and a Ajes = 0.0% sample. 

In conclusion, we have presented a novel technique 
for top-quark mass measurement in the dilepton chan- 
nel by performing in situ JES calibration using three 
different variables. Our study shows an approximately 
0.61 GeV/c 2 statistical uncertainty including overall JES 
systematic uncertainty, with 5 fb _1 LHC data. This tech- 
nique significantly improves the precision of the M top 
measurement as compared to the M top -only fit. In partic- 
ular, the overall JES uncertainty, which would be approx- 
imately 1.7 GeV/c 2 without additional calibration of JES 
with a 2% overall uncertainty assumed, is significantly 
reduced with the in situ JES calibration method. To ob- 
tain precision below 1 GeV/c 2 , one still needs to improve 
other important systematic uncertainties in the LHC ex- 
periments. However, the experience at the Tevatron pre- 
sented a well-controlled systematic uncertainty of around 
0.86 GeV/c 2 from the other systematic sources [l(|. If 
the other systematic uncertainties are controlled to the 
level of the Tevatron, we eventually can reach a precision 
of 1 GeV/c 2 in the dilepton channel. 

This work was supported by a National Research Foun- 
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FIG. 3: Plots for checking the bias in the fitted top-quark mass (a) and width of the pull distribution (b) for the 2D fit before 
corrections. With the corrections of the residual and pull width described in the text, we obtain the mass bias for various Ajes 
points at (c) for the 2D fit. 
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FIG. 4: The expected statistical uncertainty for the Mtop-only fit (a) and the Af top -AjES 2D fit (b). To evaluate the systematic 
uncertainty from the overall JES uncertainty, the mass residuals as a function of Ajes are shown for the M t0 p =173 GeV/c 2 
sample in (c) (M top -only fit) and (d) (2D fit). 
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